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Motivated by the string gas cosmological model, which predicts a blue tilt of the 
■ primordial gravitational wave spectrum, we examine the constraints imposed by cur- 

rent and planned observations on a blue tilted tensor spectrum. Starting from an 
expression for the primordial gravitational wave spectrum normalised using cosmic 
microwave background observations, pulsar timing, direct detection and nucleosyn- 



thesis bounds are examined. If we assume a tensor to scalar ratio on scales of the 
q CMB which equals the current observational upper bound, we obtain from these cur- 

' rent observations constraints on the tensor spectral index of ut < 0.79, ray < 0.53, 



and ut < 0.15 respectively. 



> 

O I- INTRODUCTION 

A stochastic background of primordial gravitational waves represents valuable informa- 
| tion about the very early universe as well as a way to discriminate between the myriad of 

O 

cosmological models currently proposed. Although a stochastic background of gravitational 



^ . waves has yet to be directly detected, efforts are being undertaken to do so by many current 
and future experiments, some of the largest being LIGO, GEO600, TAMA, and LISA 
A stochastic background of gravitational waves is characterised by the gravitational wave 
spectrum 

1_ dp gW 

where p c is the critical density of the universe, p gw is the energy density of the background 



gravitational waves and / is frequency jl]. In the above, the energy density in gravitational 
waves is written as an integral over In/, and the derivative picks out the integrand. In 
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practice, the gravitational wave spectrum is commonly assumed to depend on frequency as 
a power of /. 

Within the framework of the inflationary universe scenario the primordial gravitational 
wave spectrum is predicted to be nearly scale invariant with a slight red tilt 2|, i.e. more 
power at large scales. The reason for the red tilt is that the amplitude of the gravitational 
wave spectrum on a fixed scale k is set by the Hubble constant H at the time t^k) when 
the scale k exits the Hubble radius during the period of inflation. Smaller scales exit the 
Hubble radius later when the Hubble constant is smaller, leading to the red tilt. However, 
there do exist alternative cosmological models to the standard inflationary scenario which 
predict a blue tilt of the gravitational wave spectrum, i.e. more power at small scales, which 
have not yet been ruled out by observations. One cosmological model that predicts such a 
tilt is string gas cosmology 

String gas cosmology [3] is an approach to string cosmology which starts from the new 
degrees of freedom and symmetries which string theory contains, but particle physics-based 
models lack, namely string winding modes, string oscillatory modes and T-duality symme- 
try, and uses them to develop a new cosmological model. The claim is that by making these 
crucial stringy additions, one obtains a new cosmological model which is singularity free , 
generates nearly scale-invariant scalar metric perturbations from initial string thermody- 
namic fluctuations % |5| and provides a natural explanation for the observed dimensionality 
of space {3] (see [f| for an overview of the string gas cosmology structure formation sce- 
nario). A key result which emerges 3] from string gas cosmology is that the gravitational 
wave spectrum has a slight blue tilt, giving rise to a testable prediction different from that 
of the inflationary universe paradigm. 

Using string gas cosmology as a motivation, in this work we will focus on how a tensor 
spectrum with a blue tilt can be constrained by current observational results, and what the 
prospects for improved constraints from some planned experiments are. 

The starting point of our analysis will be an expression for the primordial gravita- 
tional wave spectrum normalised by cosmic microwave background (CMB) observations. 
Chongchitnan and Efstathiou [8j have derived just such an expression with a pivot scale 
ko = 0.002 Mpc^ 1 using the combined results from multiple surveys. Using the value 
Vs{ko) ~ 2.21 x 10~ 9 for the amplitude of the scalar power spectrum evaluated at the 
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pivot scale, they found ft gw (f) can be written as 



h 2 n gw (f) - 4.36 x l(T 15 r f£\ 



tlx 



(2) 



where fo = 3.10 x 10 18 Hz. Solving this result for the tilt of the gravitational wave spectrum, 
n T , we get the explicit expression 



In the above equations r is the tensor-to-scalar ratio evaluated at the pivot scale, 

(4) 



V T (k ) 



Vs(ko) ' 

The calculation of the tensor-to-scalar ratio depends quite sensitively on the parameters 
of the cosmological model under consideration. For that reason we choose to leave r as a 
free parameter in the main expressions calculated in this work. However, for the sake of 
examining some numerical values of the constraints derived here we will insert a value of r 
corresponding to the current upper bound into our results. Note, however, that the bounds 
we derive depend only logarithmically on r. In each case we choose to use the value of the 
tensor-to-scalar ratio given by the combined three-year Wilkinson Microwave Anisotropy 
Probe (WMAP) and lensing normalised Sloan Digital Sky Survey (SDSS) data applied to a 
typical ACDM model including tensors {2]. 

In Sections II and III we will use the expression for the tilt of the gravitational wave 
spectrum ([3]) along with several direct detection constraints to calculate the bounds on the 
required blue tilt of the spectrum. In Section IV we calculate the bound on the required blue 
tilt of the spectrum arising indirectly from the theory of big bang nucleosynthesis. In Section 
V we will investigate whether the CMB observations are compatible with the calculated 
bounds, or if they will offer even tighter constraints. We conclude with a discussion of our 
findings along with other issues related to the method used. 



II. PULSAR TIMING 



High precision measurements of millisecond pulsars provide a natural way to study low 
frequency gravitational waves. A gravitational wave passing between the earth and the 
pulsar will cause a slight change in the time of arrival of the pulse leading to a detectable 
signal. 
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A. Parkes Pulsar Timing Array 



The Parkes Pulsar Timing Array (PPTA) project [lOj is a pulsar timing experiment using 
the Parkes 64m radio telescope located in Australia with the ultimate goal of reaching the 
required sensitivity to make a direct detection of gravitational waves. The PPTA project 
hopes to make timing observations of a sample of twenty millisecond pulsars, ten or more 
of which with a precision of less than approximately 100 ns. 



Jenet et al. have developed a technique to make a definitive detection of a stochas- 
tic gravitational wave background by looking for correlations between pulsar observations. 
They have applied their method to data from seven pulsars observed by the PPTA project 
combined with an earlier data set to find a constraint on the amplitude of the characteristic 
strain spectrum. They then used this result to place a bound on the primordial gravitational 
wave spectrum 

/i 2 a, w (l/8yr) < 2.0 x 10~ 8 . (5) 

Plugging the constraint of Jenet et al. into Equation (J3j) at the frequency / = 1/8 yr ~ 
3.96 x 10~ 9 Hz we obtain a constraint 

nr < 0.0477tof 4 - 59xl0 'V (6) 



The WMAP+SDSS data places a bound r < 0.30 on the tensor-to-scalar ratio [2|. Insert- 
ing r = 0.30 into the above equation we find that the current pulsar timing observations 
constrain the blue tilt of the tensor spectrum to ny < 0.79. 

Jenet et al. have also used simulated data to determine the upper bound on the primordial 
gravitational wave spectrum expected from future pulsar observations. Using a simulated 
data-set of twenty pulsars timed with an rms timing residual of 100 ns over 5 years they 



calculated 



12| 



h 2 Q gw (l/8jT) < 9.1 x HT 11 . (7) 



Plugging this improved constraint into Equation ([3]) at the frequency / = 1/8 yr we get a 
bound 

„ T < 0.0477* (8) 

and again using the value r = 0.30, we find that, in the absence of a detection, future pulsar 
timing observations could tighten the constraint on the blue tilt to nx < 0.53. 
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III. INTERFEROMETERS 



Interferometer experiments offer a way to directly measure the gravitational wave strain 
spectrum with many observatories currently running or planned for the future. Interferom- 
eters in different locations form a network that will search for a correlated signal between 
detectors beneath uncorrelated detector noise, in order to improve sensitivity. 



A. LIGO 

The Laser Interferometer Gravitational Wave Observatory (LIGO) |l3j is a ground based 
interferometer project operating in the frequency range of 10 Hz - a few kHz. LIGO consists 
of two collocated Michelson interferometers in Hanford, Washington, HI with 4 km long 
arms, and H2 with 2 km long arms, along with a third interferometer in Livingston Parish, 
Louisiana, LI with 4 km long arms. 

Most recently LIGO has performed its fourth science run, S4, with improved interferome- 
ter sensitivity. Abbott et al. [k^] have used the S4 data to calculate a limit on the amplitude 
of a frequency independent gravitational wave spectrum. They found a bound 

tt gw < 6.5 x 1(T 5 (9) 

in the frequency range 51-150 Hz. Inserting this value into Equation ([3]) at the frequency 
/ = 100 Hz we get a constraint 

/1.49 x 10 10 h 2 \ , . 

n T < 0.0223 In ( J . (10) 



The WMAP+SDSS data also provides a value of h = 0.716 for the Hubble parameter 
Inserting this along with r = 0.30 into the above bound, we find that the current LIGO 
results place a constraint n-p < 0.53 on the blue tilt of the tensor spectrum. 

The final phase of LIGO, named Advanced LIGO, hopes to reach a detection sensitivity 



of 



13| 



fi^-io- 9 . (11) 

Plugging this value into Equation (j3J) at / = 100 Hz, we find that if Advanced LIGO does 
not make a positive detection of a gravitational wave background then it will place a bound 
on the blue tilt of the tensor spectrum 

/2.29 x 10 5 /i 2 \ , s 

n T < 0.0223 In . (12) 
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Using r = 0.30 and h = 0.716 in this expression, Advanced LIGO would then constrain the 
blue tilt of the tensor spectrum to n T < 0.29. 



B. LISA 

The Laser Interferometer Space Antenna (LISA) [2| is a planned space-based interferom- 
eter experiment operating in the mHz range. LISA will consist of three drag-free spacecraft 
each at the corner of an equilateral triangle with sides of length 5 x 10 9 m. Each spacecraft 
has two optical assemblies pointed towards the other two spacecraft forming three Michelson 
interferometers. This triangle formation will orbit the sun in an Earth-like orbit separated 
from us by approximately fifty million kilometers. The goal of LISA is to reach a sensitivity 



of 



15| 

h 2 tt gw (lmRz) ~ 1 x 10~ 12 . (13) 



At the LISA sensitivity level one would expect gravitational wave signals from super-massive 
black hole binaries, other binary systems and super-massive black hole formation to be 
present. Assuming these predicted signals could somehow be removed and LISA does not 
detect any primordial signal, we can plug this predicted limit into Equation dSJ) at / = lmHz 
to obtain a limit on the blue tilt of the primordial gravitational wave spectrum 

n T < 0.0299 In ( 229 x 10 " \ . ( 14) 



Inserting r = 0.30 into this equation we find that LISA could potentially place a constraint 
n T < 0.20 on the blue tilt. 



IV. NUCLEOSYNTHESIS 



The theory of big-bang nucleosynthesis (BBN) successfully predicts the observed abun- 
dances of several light elements in the universe. In doing so, BBN places constraints on a 
number of cosmological parameters. This in turn results in an indirect constraint on the 
energy density in a gravitational wave background as follows: the presence of a significant 
amount of gravitational radiation at the time of nucleosynthesis will change the total en- 
ergy density of the universe, which affects the rate of expansion in that era, leading to an 
over-abundance of helium and thus spoiling the predictions of BBN [l|. Assuming N u = 4.4, 
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where N v is the effective number of neutrino species at the time of nucleosynthesis, the BBN 



bound is 



13| 



[ h n gw (f)d(\nf) < 1.5 xl(T 5 . (15) 
Jh 

Plugging Equation (j2J) into the left-hand side and performing the integration we obtain the 
inequality 

h ^<3.4xl0 9 ^^. (16) 

tit r 

In order to apply the above result, we must discuss the two integration limits fx and f 2 - 
The lower cutoff frequency fi corresponds to the Hubble radius at the time of BBN and takes 
the value j\ ~ 10~ 10 Hz. For wavelengths larger than the Hubble radius, the gravitational 
waves are frozen out [l6] (see e.g. jl^] for a review) and thus do not act like radiation. The 
upper cutoff frequency f 2 is the ultraviolet cutoff. We will take it to be given by the Planck 
frequency, i.e. f 2 = fpi = 1.86 x 10 43 Hz. Substituting these two limits into Equation (TIB"]) 
along with the WMAP+SDSS values r = 0.30 and h = 0.716 then solving numerically for 
ny, we find the bound on the blue tilt of the tensor spectrum from BBN to be 

n T <0.15. (17) 

Had we instead inserted for f 2 the scale of grand unification, 10 16 GeV, or the Hubble rate 
during a simple large field inflation model, which is 10 13 GeV, the bound would be slightly 
relaxed to 0.16 or 0.17, respectively. Thus, the dependence of the bound on the uncertain 
ultraviolet cutoff scale f 2 is quite mild. 

Lastly, we do not want Q gw > 1 at any scale within the integration bounds. Since we are 
working with such large frequencies we should check to be sure that this condition is satisfied. 
Substituting the value of the tensor spectral index determined by BBN back into Equation 
(J2]) we find that fl gw (fpi) = 3.93 x 10~ 6 , meaning our requirement is indeed satisfied for all 
frequencies within the interval of integration. 



V. COSMIC MICROWAVE BACKGROUND 

Observations of the CMB have implications for a wide variety of topics including con- 
straining inflation, dark matter and large-scale structure. The mission of WMAP is to 
produce full-sky maps of the CMB anisotropy and the recently published three-year results 
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are an improvement upon previous observations. A reduction in instrument noise produced 
spectra which are three times more sensitive in the noise limited region, independent years of 
data allow for cross-checks, the instrument calibration and response have been better char- 
acterised and a thorough analysis of the polarisation data has improved the understanding 
of the data [18 



Using three year WMAP data, the derived angular power spectrum of the 
temperature anisotropy, Cf T , where I is the multipole moment, is cosmic variance limited 
to I = 400 and the signal to noise ratio exceeds unity to / = 1000 This high precision 
cosmological data provides another method of placing constraints on the value of the tensor 
spectral index. 

Using the Code for Anisotropies in the Microwave Background (CAMB) [l9] we can 
simulate how a blue tilt of the primordial gravitational wave background would effect the 
anisotropies in the CMB. To examine possible constraints we employ the following method: 
first, we calculate Cf T using CAMB for each of the three current bounds on calculated 
in the previous sections; second, we calculate Cf T using CAMB, this time with a standard 
inflationary relation: tlt = — r/8 2]; finally, we compare the output Cj T data for the models 
with a blue tilt against the output for the model with the "usual" value of the tensor spectral 
index (the above relation from inflationary cosmology). For consistency with the previous 
sections, when running CAMB we choose our input cosmological parameters to be those 
calculated using the WMAP+SDSS data for a ACDM model with tensors jsj]. 

From Figured] we can clearly see that the power spectrum of the temperature anisotropy 
for models with a blue tensor spectral index does not vary much from that calculated using a 
standard inflationary definition of the tensor spectral index. In fact, the difference is within 
the cosmic variance error at all I < 1000 for each of the three bounds calculated using the 
PPTA observations, LIGO observations and the theory of BBN. Thus, we find that the 
cosmic microwave background does not offer any tighter constraints on the blue tilt of the 
gravitational wave spectrum than those already calculated. 



VI. CONCLUSIONS 



A stochastic background of primordial gravitational waves is a prediction of many cosmo- 
logical models. Assuming the gravitational wave spectrum depends as a power on frequency, 
then this spectrum can be characterised by its tilt and amplitude. Most models predict the 
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FIG. 1: Magnitude of the difference between the the TT power spectrum for a model with a blue 
tensor spectral index and the TT power spectrum for a model with a tensor spectral index defined 
as: ny = — r/8. Shown here are the cases bt = 0.15 (green), tit = 0.53 (yellow) and ut = 0.79 
(orange) respectively. The dashed line represents the cosmic variance error at each I. 



tilt to be nearly scale-invariant but slightly red, while some models, like string gas cosmology, 
predict a slight blue tilt. Although a gravitational wave background has yet to be directly 
detected, observational results can already be used to constrain it. 

Using the current results from pulsar timing observations, direct detection observations 
and the theory of nucleosynthesis, we have placed bounds on the possible blue tilt of the 



20j 



m 



gravitational wave background. After completion of this work we became aware of 
which a master equation was derived which relates the short wavelength observable £l gw (f) 
to the tensor to scalar ratio measured with the CMB. The goal of that work was to develop 
a formulation which is as general as possible. In particular, the equation of state parameter 
and the tensor spectral index are taken to be arbitrary functions of the scale factor and 
wavenumber respectively, not constants as is often assumed. Their master equation is thus 
a more general version of our "master" equation (|2J) and we could have just as easily used 
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it as our starting point. In fact, we have confirmed that by choosing a constant value 
w = 1/3 for the equation of state parameter (which is described as the most logical value in 
[20I ]) and numerical values for other cosmological parameters that match our choices above, 
we can indeed re-derive all of our results from the master equation of [20j]. We consider 
this a good consistency check for the constraints derived in this work. The authors of [20J 
do include a discussion of some of the same types of observations mentioned here, namely 
laser interferometer and pulsar timing, however, we stress that they do not use the current 
numerical constraint from any particular observatories to compute actual upper bounds on 
Ut, which was the purpose of this work. The authors of 2Q] also discuss a constraint on tit 
coming from BBN, but they take the constraint on Q gw from BBN to be a constant across 
all frequencies rather than integrating their master equation as was done in this work. In 
the end they find a weaker bound on the tilt, < 0.36, than the one obtained here. 

By far the tightest constraint on the tilt comes from big-bang nucleosynthesis. If we take 
the tensor to scalar ratio on CMB scales to be given by the current observational upper 
bound, and if we take the ultraviolet cutoff scale in the spectrum of gravitational radiation 
to be the Planck scale, then the bound is tit < 0.15, tighter than even Advanced LIGO, 
the future PPTA, and LISA can hope to achieve. It is not surprising that nucleosynthesis 
provides the tightest bounds on a blue tilt of the gravitational wave spectrum since nucle- 
osynthesis probes physics on scales much smaller than the other experiments we analysed, 
and spectra with blue tilts have more power on the smallest scales. That is, BBN gives us 
the largest "lever arm" to probe gravitational waves in conjunction with CMB observations, 
a point also made in 21J. From our results we can clearly see the trend that the bound on 
the blue tilt of the tensor spectral index tightens as the length scale probed by the given 
experiment decreases. 

Simulations of the angular power spectrum of the temperature anisotropics in the CMB 
did not offer any tighter constraints on the tensor spectral index, with each of the constraints 
calculated in this paper producing a temperature power spectrum that was within the cosmic 
variance error of one calculated for a standard ACDM model with tensor modes included. 

As mentioned in the Introduction, Equation (j2J) has been normalised at the scale of 
cosmic microwave background observations. However, those experiments probe scales that 
are approximately ten orders of magnitude larger than those probed by the PPTA and 
approximately nineteen orders of magnitude larger than those probed by LIGO, with LISA 
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probing between the two. Extrapolating between such a large difference in scales is not 
straightforward and we should note that in [8| the authors conclude from their analysis 
that even within the framework of the inflationary universe paradigm the formula for the 
primordial gravitational wave spectrum (j2J) is too restrictive, and they believe it is indeed 
not possible to extrapolate reliably over such a large difference in scales. Whether or not this 
is the case in the string gas cosmology model should perhaps be examined more carefully in 
future work. 

Continuing with string gas cosmology, we conclude that the current bounds on the tilt of 
the gravitational wave spectrum are weak. The predicted magnitude of the blue tilt of the 
gravity wave spectrum is thought to be comparable to the magnitude of the red tilt of the 
spectrum of scalar metric fluctuations [7]. If the latter is taken to agree with the current 
bounds, we predict a blue tilt of less than = 0.1 which will not be easy to detect. There 
may, however, be models similar to string gas cosmology in which the scalar and tensor tilts 
are not related, and for which planned experiments could set valuable constraints on the 
model parameter space. 
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